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ABSTRACT 

Following the recent abundance measurements of Mg, Al, Ca, Fe, and Ni in the black 
hole X-ray binary XTE J1118-I-480 using medium-resolution Keck II/ESI spectra of the 
secondary star (Gonzalez Hernandez et al. 2006), we perform a detailed abundance 
analysis including the abundances of Si and Ti. These element abundances, higher than 
solar, indicate that the black hole in this system formed in a supernova event, whose 
nucleosynthetic products could pollute the atmosphere of the secondary star, providing 
clues on the possible formation region of the system, either Galactic halo, thick disk, 
or thin disk. We explore a grid of explosion models with different He core masses, 
metallicities, and geometries. Metal-poor models associated with a formation scenario 
in the Galactic halo provide unacceptable fits to the observed abundances, allowing us 
to reject a halo origin for this X-ray binary. The thick-disk scenario produces better 
fits, although they require substantial fallback and very efficient mixing processes be- 
tween the inner layers of the explosion and the ejecta, making quite unlikely an origin 
in the thick disk. The best agreement between the model predictions and the observed 
abundances is obtained for metal-rich progenitor models. In particular, non-spherically 
symmetric models are able to explain, without strong assumptions of extensive fallback 
and mixing, the observed abundances. Moreover, asymmetric mass ejection in a super- 
nova explosion could account for the required impulse necessary to launch the system 
from its formation region in the Galactic thin disk to its current halo orbit. 

Subject headings: black holes: physics — stars: abundances — stars: evolution — stars: indi- 
vidual (XTE Jll 18+480) — supernovae: general — X-rays: binaries 
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1. Introduction 

The low-mass X-ray binary XTE Jl 118-1-480 
is the first identified black hole moving in 
Galactic halo regions (Wagner et al. 2001; 
Mirabel et al. 2001). Since it was discovered 
during a faint outburst on UT 2000 March 
29 (Remillard et al. 2000), it has been inten- 
sively studied in both the X-ray and optical 
spectral regions. During the decay of the out- 
burst, McClintock et al. (2001) and Wagner 
et al. (2001) determined the radial-velocity 
curve of the companion star, yielding a mass 
function /(M) ~ 6 M.q. The companion star 
was classified as a late-type main-sequence 
star with a mass of 0.1-0.5 M© (Wagner et al. 
2001). 

By modelling the light curve, McClintock 
et al. (2001) derived a lower limit to the or- 
bital inclination, i > 55°, and consequently 
an upper limit to the black hole mass of 
-^BH ^ 10 Mq. Additional evidence for a 
high inclination (i > 60°) comes from mea- 
surements of tidal distortion (Prontera et al. 
2001), whereas the lack of dips or eclipses for a 
Rochc-lobc filling secondary yields upper lim- 
its of i > 80° and Mrh > 7.1 Mq. Later, 
Gelino et al. (2006) derived an orbital inclina- 
tion of 68° ± 2° , by modeling the optical and 
infrared ellipsoidal light curves of the system 
in quiescence. This value of the inclination al- 
lowed them to better constrain the black hole 
mass at Mbh = 8.53 ± 0.60 M©. 

The system is placed in the Galactic halo, 
with an extraordinarily high Galactic latitude 
(6 62.3°), and a height of ~1.6 kpc above 
the Galactic plane, according to its distance of 
1.85±0.36 kpc (Wagner et al. 2001). This ap- 
pears surprising since all other black hole bi- 
naries are located in the Galactic disk. An ac- 
curate measurement of its proper motion cou- 
pled with its distance provides space-velocity 
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components {U, V) which seem consistent 
with those of some old halo globular clus- 
ters (Mirabel et al. 2001). This opened the 
possibility that the system originated in the 
Galactic halo, and therefore, that the black 
hole could be either the remnant of a super- 
nova (SN) in the very early Galaxy or the 
result of direct collapse of an ancient mas- 
sive star. However, the galactoccntric orbit 
crossed the Galactic plane many times in the 
past, and the system could have formed in 
the Galactic disk and been launched into its 
present orbit as a consequence of the "kick" 
imparted during the SN explosion of a massive 
star (Gualandris et al. 2005) . Recent observa- 
tions with the 10-m Keck II telescope revealed 
that the secondary star has a supcrsolar sur- 
face metalhcity ([Fc/H] = 0.2 ± 0.2, Gonzalez 
Hernandez et al. 2006), confirming the ori- 
gin of the black hole in a SN event. Thus, 
if the system originated in the Galactic halo, 
the element abundances of the secondary star 
must have been enriched by a factor of 5-25 
depending on whether its initial metallicity re- 
sembled a thick-disk star or a halo star. 

Element abundances of secondary stars of 
X-ray binaries have been studied for the sys- 
tems Nova Scorpii 1994 (Israelian et al. 1999; 
Gonzalez Hernandez et al. 2007), A0620-00 
(Gonzalez Hernandez et al. 2004), Centau- 
rus X-4 (Gonzalez Hernandez et al. 2005), 
XTE J1118-F480 (Gonzalez Hernandez et al. 
2006, hereafter Paper I), and V4641 Sagittarii 
(Orosz et al. 2001; Sadakane et al. 2006). Ah 
of these X-ray binaries show metallicities close 
to solar independent of their location with re- 
spect to the Galactic plane, and possible sce- 
narios of pollution from a SN or hypernova 
have been discussed. In this paper, we com- 
pare in detail different scenarios of the possi- 
ble enrichment of the secondary star from SN 
yields, providing conclusions on the formation 
region (Galactic halo, thick disk, or thin disk) 
of this halo black hole X-ray binary. 
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Fig. 1. — Radial velocities of XTE J1118+480 folded 
on the orbital solution of the data, together with the 
best-fitting sinusoid. Individual velocity uncertainties 
are < 7 km s~^ and are not plotted because they are al- 
ways smaller than the symbol size. The bottom panel 
shows the residuals of the fit. 

2. Observations 

As already reported in Paper I, we ob- 
tained 74 medium-resolution spectra of the 
black hole X-ray binary XTE J1118-F480, in 
quiescence, with the Echelle Spectrograph and 
Imager (ESI; Sheinis et al. 2002) at the 10-m 
Keck II telescope on UT 14 February 2004. 
The data covered the spectral range 4000- 
9000 A at a resolving power A/AA « 6,000. 
We also observed ten template stars with 
spectral types in the range K0V--M2V with 
the same instrument and spectral configura- 
tion. The exposure time was fixed at 300 s to 
minimize the eff^ects of orbital smearing which, 
for the orbital parameters of XTE J1118-I-480, 
is in the range 0.6-26.6 kms~^, smaller than 
the instrumental resolution of 50 kms~^. All 
of the spectra were reduced in a standard 
manner. 

3. Revised Orbital Parameters 

We extracted the radial velocities by cross- 
correlating each target spectrum with the 



spectrum of a K5V template star, using the 
software MOLLY developed by T. R. Marsh. 
We fitted these data with a sine wave using a 
x2 method providing the following orbital so- 
lution (see Fig.d]): 7 = 2.7±1.1 kms"\ K2 = 
708.8 ± L4 km s-\ P = 0.16995 ± 0.00012 d, 
and To = 2, 453, 049.93346 ± 0.00007 d, where 
To is defined as the corresponding time of 
the closest inferior conjunction of the com- 
panion star, and the quoted uncertainties are 
la. This orbital period, P, and the velocity 
amplitude of the orbital motion of the sec- 
ondary star, K2, lead to a mass function of 
/(M) = 6.27 ± 0.04 M0, consistent with (but 
more precise than) previous results (McClin- 
tock et al. 2001; Wagner et al. 2001; Torres 
et al. 2004). 

The derived radial velocity of the center of 
mass of the system agrees somewhat (at the 
3a level) with previous studies (7 = -|-26 =b 17 
kms~\ McClintock et al. 2001; 7 = -15 ±10 
kms~^, Wagner et al. 2001; 7 = +16 ± 6 
kms~^, Torres et al. 2004). Note that our 
medium-resolution data have a factor of > 4 
higher spectral resolution than the data these 
authors used. We estimate the uncertainty in 
our individual radial-velocity measurements 
of typically 6 kms~^ (see Fig. [1]). 

3.1. Secondary Spectrum 

The individual spectra were corrected for 
their radial velocity and combined in order to 
improve the signal-to-noise ratio (S/N). After 
binning in wavelength in steps of 0.3 A, the 
final spectrum had S/N ~ 80 in the contin- 
uum in the red spectral region. This spec- 
trum, displayed in Fig. [21 was compared with 
ten template stars having spectral types KOV 
to M2V. The best fit shows a K5V star rather 
than the later spectral types (K7/8V or even 
M) suggested in previous studies (Wagner et 
al. 2001; Torres et al. 2004). 

We should remark that the spectrum of the 
secondary star in this system and the spectra 



3 



«5 Hy 



Hell? 




XTEJl 118+48 



Hel 



Hel Hel 



Can 



f Ha Lil 



6000 7000 
>. (ANGSTROMS) 



Template star: BD-053763 (K5V) 




Lil 



6000 7000 
X (ANGSTROMS) 



Template star: BD+52857 (K7/8V) 



*,iJ 



l (ANGSTROMS) 



Fig. 2. — Observed spectrum of the seeondary star of XTE Jl 118+480 (top panel) and of two properly broadened 
templates (BD— 053763, middle panel; BD+52857, bottom panel). 



of the template stars were normalized using 
the same procedtirc. The continuum was fit- 
ted with a low-order spline, in order to avoid 
the smoothing of possible existing broad TiO 
bands. Following Marsh et al. (1994), we 
computed the optimal value of v sin i by sub- 
tracting broadened versions of the K5V tem- 
plate (in steps of 1 kms~^) and minimizing 
the residual. We used a spherical rotational 
profile with linearized limb-darkening e = 0.8 
(Al-Naimiy 1978) due to the spectral type K 
of the secondary star. The best fit corre- 
sponds to V sin i = lOOlf^ kms~^, where the 
uncertainties have been derived by assuming 
extreme cases for e = — 1. Our derived ro- 
tational velocity, combined with our value of 
the velocity amplitude, K2, implies a binary 



mass ratio q = 0.027 ± 0.009, in agreement 
with previous resTilts (Torres et al. 2004, and 
references therein). 

4. Chemical Analysis 
4.1. Stellar Parameters 

The normalized spectra of X-ray transients, 
although observed in quiescence, show appar- 
ently weaker stellar lines of the secondary star 
due to the veiling caused by the accretion disk. 
The veiling from the accretion disk was esti- 
mated to be ~ 65% ±8% in the spectral range 
5800-6400 A in December 2000 and January 
2001, and ~ 40% ± 10% in January 2003, 
by performing standard optimal subtraction 
techniques with K5V-M0V template stars. 
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Table 1 

Ranges and Steps of Model Parameters 



i'aiaui(n(*r 


Range 




Tcfi 


3500 5000 K 


100 K 


log [9/ (cm s2)] 


4^5 


0.1 


[Fe/H] 


-1.5 1 


0.1 


/4500 


0^2 


0.1 


mo 


-0.00091 


-0.00010 



I 0.6: 

O 0.4- 

o 



n.n ro 

4000 4200 



4400 4600 



4800 5000 




Fig. 3. — Distributions obtained for each parameter using Monte Carlo simulations. The labels at the top of each 
bin indicate the number of simulations consistent with the bin value. The total number of simulations was 1000. 



As shown in Paper I, we tried to infer 
the stellar parameters, Tcfj and log 5, and the 
metallicity [Fe/H], of the secondary star tak- 
ing into account the veiling from the accre- 
tion disk, defined as a linear function of wave- 



length and thus described with two additional 
parameters, the veiling at 4500 A, /4500 = 
-^dfsk7^stc°°> and the slope, mo- Note that the 
total flux is defined as Ftotai = -^disk + 
where Fdisk and Fgec are the flux contribu- 
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Fig. 4. — Best synthetic spectral fits to tlie Keck/ESI spectrum of tfie secondary star in tfie XTE J1118+480 system 
(bottom panel) and the same for a template star (properly broadened) shown for comparison (top panel). Synthetic 
spectra are computed for solar abundances (dashed line) and best-fit abundances (solid line). 
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Fig. 5. — The same as in Fig. 4, but for the spectral range 5510-5550 A. 
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Fig. 6. — The same as in Fig. 4, but for the spectral range 6420-6480 A. 



tions of the disk and the continuum of the sec- 
ondary star, respectively. This procedure in- 
volves a x2 minimization routine which com- 
pares several features of the stellar spectrum 
with a grid of synthetic spectra computed us- 
ing the LTE code MOOG (Sneden 1973). We 
used a grid of LTE model atmospheres (Ku- 
rucz 1993) and atomic line data from the Vi- 
enna Atomic Line Database (VALD, Piskunov 
1995). The oscillator strengths of relevant 
lines were adjusted until they reproduced the 
solar atlas (Kurucz et al. 1984) with solar 
abundances (Grevesse et al. 1996). The 
changes we applied to the log gf values taken 
from the VALD database were Alogg/ < 0.2 
dex. 

We selected nine spectral features contain- 
ing in total 30 lines of Fe I and 8 lines of Ca I 
with excitation potentials between 1 and 5 eV. 
The five free parameters were varied in the 
ranges given in Table [H For each given iron 
abundance in the range [Fe/H] < 0, the Ca 



abundance was fixed according to the Galac- 
tic trend of Ca (Bensby et al. 2005), while for 
[Fe/H] > 0, we assumed [Ca/Fe] = 0. A rota- 
tional broadening of 100 kms~^ and a limb 
darkening e = 0.8 were adopted. The micro- 
turbulence, ^, was computed using an exper- 
imental expression as a function of effective 
temperature and surface gravity (Allende Pri- 
eto et al. 2004). 

The result, already presented in Paper I, 
provides as most likely values T^g = 4700 ± 
100 K, log[5/cm s^] = 4.6 ± 0.3, [Fe/H] = 
0.18 ± 0.17, /4500 = 0.85 ± 0.20, and mo = 
-0.0002 ± 0.0001. The la uncertainties of 
the five free parameters were determined us- 
ing 1000 realizations whose corresponding his- 
tograms are displayed in Fig. [3l 

The stellar parameters derived, especially 
the effective temperature, could provide im- 
portant implications on the Gelino et al. 
(2006) determination of the orbital inclination 
of the system (i ~ 68°). These authors mod- 
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Table 2 

Abundance Uncertainties in XTE J1118+480 



Element 


[X/H]lte" 








Avoil 




^lincs 


Mg 


0.35 


0.12 





-0.10 


0.20 


0.25 


1 


Al 


0.60 


0.12 


0.05 





0.15 


0.20 


1 


Si 


0.37 


0.03 


-0.07 


0.07 


0.18 


0.21 


2 


Ca 


0.15 


0.03 


0.13 


-0.16 


0.11 


0.23 


5 


Ti 


0.32 


0.18 


0.12 


-0.03 


0.14 


0.26 


3 


Fe 


0.18 


0.08 


0.06 


0.04 


0.13 


0.17 


5 


Ni 


0.30 





0.10 


0.12 


0.14 


0.21 


2 


Lid 


1.78 


0.12 


0.15 


0.05 


0.15 


0.25 


1 



^Element abundances of tiie secondary star (calculated assuming LTE) are 
[X/H] = log[Ar(X)/Af (H)]star -log[Af(X)/Af (H)]sun, where Af (X) is the number 
density of atoms. Uncertainties, A[X/H], are at the la level and take into 
account the uncertainties in the stellar and veiling parameters. 

''The total error was estimated as Atot = ,/A2 + + A? + A^ ... 

"jui. Y ^ Toff log g veil 

'^Number of features analyzed for each element. 

''Li abundance is expressed as loge(Li)NLTE = log[A'^(Li)/A''(H)]iMLTE + 12. 



Note. — The uncertainties from the dispersion of the best fits to different 
features, A^, are estimated using the following formula: A^ = <y/\/N, where 
cr is the standard deviation of the measurements. 



eled the optical and infrared (IR) ellipsoidal 
light curves of XTE J1118+480 in quiescence, 
assuming a K7V spectral type {T^s ~ 4250 K) 
for the secondary star. However, our spectro- 
scopic value (Teff 4700 K) would require 
more contribution of the flux from the accre- 
tion disk in the K band. Gelino et al. (2001) 
derived an inclination of z ~ 41°, yielding 
a black hole mass of Mbh ~ H Mq for the 
A0620-00 system, by adopting Teg 4600 K, 
which is 300 K lower than the spectroscopic 
value reported by Gonzalez Hernandez et al. 
(2004). Hynes et al. (2005) suggested that 
this different effective temperature of the 
secondary star in A0620-00 would require 
a larger disk contribution in the K band, 
and therefore a higher inclination. In fact, 
Gelino et al. (2001) commented that if the 
i^-band disk veiling as high as 50% of the 
total flux, the derived inclination would in- 
crease to 60°. The black hole mass would 



then drop to Mbh ~ 5 M©. Although this 
is an extreme case, milder IR veiling could 
still have a substantial impact on the derived 
black hole mass. Similarly, the black hole 
mass in XTE J1118-I-480 might be signifi- 
cantly affected when using our value of the 
effective temperature. 

4.2. Stellar Abundances 

We inspected several spectral regions in the 
observed Keck/ESI spectrum of the secondary 
star, searching for suitable lines for a detailed 
chemical analysis. Using the derived stellar 
parameters, we first determined the Fe abun- 
dance by comparing synthetic spectra with 
each individual feature in the ESI spectrum 
(see Table E]). In Fig. g] (here) and Fig. 1 of 
Paper I, we display some of the spectral re- 
gions analyzed to obtain the Fe abundance. 
This figure also shows the best synthetic spec- 
tral fit to the observed spectrum of a template 
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Fig. 7. — Abundance ratios of the secondary star in XTE J1118+480 (blue wide cross) in comparison with the 
abundances of G and K metal-rich dwarf stars. Galactic trends were taken from Gilli et al. (2006). The size of the 
cross indicates the uncertainty. Filled and empty circles correspond to abundances for exoplanet host stars and stars 
without known exoplanet companions, respectively. The dashed-dotted lines indicate solar abundance values. 



Star (HIP 17420 with Tcs= 4801 K, log g= 
4.633, and [Fe/H]= -0.14 dex) using the stel- 
lar parameters and abundances determined 
by Allende-Prieto et al. (2004). We only 
use as abundance indicators those features 
which were well reproduced in the template 
star. The chemical analysis is summarized 
in Table [2j The errors in the element abun- 
dances show their sensitivity to the uncertain- 
ties in the effective temperature {At^^), grav- 
ity (Alog^), veiling (Aveiiing)) and the disper- 
sion of the measurements from different spec- 
tral features (Aq-). In Table[2]we also state the 
number of features analyzed for each element. 

The abundances of Ti and Si were mainly 
derived from several lines in the spectral re- 
gion 5920-5960 A where some telluric lines 
are present. However, since 74 spectra with 



different radial velocities in the range ±710 
km s^^ were combined to generate the average 
spectrum of the secondary star, these telluric 
lines must have been smoothed out. 

The Mg abundance was derived from one 
spectral line (see Fig. [5]) and the error as- 
sociated with the dispersion of the measure- 
ments, cj, was assumed to be the average dis- 
persion of Fe, Ca, and Ni abundances, and 
in this case, A^- = a. The same prescription 
was adopted for the analysis of Al and Li (see 
Fig. 1 of Paper I). The best fit to the Li I 
6708 A feature provides an LTE abundance 
of loge(Li)LTE = 1-61 ± 0.25. We estimated 
the non-LTE abundance correctioE0 for this 
element from the theoretical LTE and non- 
LTE curves of growth in Pavlenko & Magazzii 

"^AnLTE = loge(X)NLTE — loge(X)LTE. 
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Table 3 

Element Abundance Ratios in XTE Jl 118+480 



Element [X/FelxxElllS ^[X/Fe],XTE1118 [^/-^'^'^l stars '^stars A,,,stars 



Mg 


0.17 


0.22 


0.02 


0.08 


0.02 


Al 


0.42 


0.15 


0.14 


0.10 


0.02 


Si 


0.19 


0.17 


-0.03 


0.04 


0.01 


Ca 


-0.03 


0.23 


-0.13 


0.05 


0.01 


Ti 


0.14 


0.22 


0.05 


0.08 


0.01 


Ni 


0.12 


0.12 


-0.02 


0.05 


0.01 



*Uncertainties in the element abundance ratios ([X/Fe]) in the secondary star in 
XTE J1118+480. 



Note. — [^/Fc]j,jj^j.j, indicate the average values of 24 stars with iron content in the 
range 0.01 to 0.35 corresponding to Ict in the [Fe/H] abundance of the secondary 
star in XTE J1118+480, taken from Gilli et al. (2006). The uncertainty in the 
average value of abundance ratios in the comparison sample is obtained as A^^gtars = 
"■stars/x/JV, where (Tstars is the standard deviation of the measurements and N is the 
number of stars. 



Table 4 

Supernova/Hypernova Model Parameters 



Mms 




E51 


Z 


M(Fe) 


Mcut 






Abundance Pattern 


40 


15.8 


1 


0.001 


5.8982E-01 


2.51 


5.56 


8.07 


Fig. 9a 


40 


15.8 


1 


0.001 


2.6581E-01 


3.01 


5.06 


8.07 


Fig 


. 9a 


40 


15.8 


30 


0.001 


7.6394E-01 


4.05 


4.02 


8.07 


Fig, 


, 9b 


40 


15.8 


30 


0.001 


9.2867E-02 


5.03 


3.04 


8.07 


Fig. 


, 9b 


40 


15.8 


1 


0.004 


3.7283E-01 


2.50 


5.57 


8.07 


Fig. 


10a 


40 


15.8 


1 


0.004 


2.6649E-02 


3.03 


5.04 


8.07 


Fig. 


10a 


40 


15.8 


30 


0.004 


9.8590E-01 


3.03 


5.04 


8.07 


Fig. 


10b 


40 


15.8 


30 


0.004 


2.0354E-01 


4.02 


4.05 


8.07 


Fig. 


10b 


30 


11 


1 


0.004 


1.4667E-01 


2.52 


5.49 


8.01 


Fig. 


11a 


30 


11 


1 


0.004 


1.7806E-03 


3.04 


4.97 


8.01 


Fig. 


11a 


30 


11 


20 


0.004 


5.7203E-01 


3.04 


4.97 


8.01 


Fig. 


lib 
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Fig. 


13b 



''Mbh i is the mass of the remnant after the explosion and before the secondary started to 
transfer matter onto the compact object. Mgjj j = 8.53 ± 0.60 Mq is the observed mass of the 
black hole. 

Note. — Supernova and hypernova explosion models used in Figs. 9-13. The quantities shown 
arc the main-sequence mass, the mass of the He core, the explosion energy, £51 = Ek /^O"^ erg, 
the ejected Fe mass, the metallicity of the model, the mass-cut, the mass of the fallback matter, 
and the final remnant mass. The mass are in units of Mq. 
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(1996). We found Anlte = 0.17. Due to the 
weakness of the absorption we consider this 
abundance estimate given in Table [2] as an 
upper limit. 

5. Discussion 

As aheady discussed in Paper I, the Fe 
abundance of the secondary star is slightly 
higher than solar, but similar to that of many 
stars in the solar neighborhood. The abun- 
dances of other elements listed in Table [2] rel- 
ative to iron are compared in Fig. [7] with the 
Galactic trends of these elements in the rele- 
vant range of metallicities. Moderate anoma- 
lies are found only for Al. In Table [3] we show 
the element abundance ratios in the secondary 
star in XTE J1118+480 and the average val- 
ues in stars with iron content in the range 
0.01 < [Fe/H] < 0.35, the comparison sample, 
corresponding to a la uncertainty in the iron 
abundance of the companion star. Whereas 
Ca and Ti are consistent with the average 
values of the comparison sample, Ni and Si, 
at la, and especially Al, at 2a, appear to be 
more abundant than the average values of the 
stars in the comparison sample. 

The present location and space- velocity 
components {U, V) of the system might sug- 
gest that the system belongs to the Galactic 
halo, but the derived metallicity makes this 
possibility less likely. One could include the 
metallicity in the expression given in Bensby 
et al. (2005) to estimate the relative like- 
lihoods that a star belongs to the Galactic 
thin disk, thick disk, and halo. The equations 
could be written as follows: 



P -f ^ 

J^thin— disk JD p ■) 



Pi^ KiX exp 



Pf 



thick— disk 



TD- 



-Ptd 

p '' 



halo 



/h 



asym,? y 



X exp 



2ct2(7. la^Vi 
M/2LSR ([Fe/H] - [Fe/H]^,^^ J2 ' 
2a2wi 2cr2[Fo/H]i 

P - /dPd + /tdPtd + /hPh 

and Kj 



P 



(27r)2(T;7.crvACrvi/,cr[Fo/H], ' 

where the subscript i indicates the three pop- 
ulations D (thin disk), TD (thick disk), and 
H (halo). The total probability, P, takes into 
account the fraction of stars belonging to each 
population in the solar neighborhood (/d = 
0.94, /td = 0.06, and /h = 0.0015; Bensby 
et al. 2003). The velocity distributions of 
each population with respect to the local stan- 
dard of rest (LSR) are centered at zero except 
for the component Flsr, whose center is dis- 
placed according to Ksym.j (with Vasym,D = 

— 15kms~^, Vasym.TD = —46 kms~^, and 
Kisym,_ff = —220 kms~^). The metallic- 
ity distributions have been characterized 

with [Fe/H],,y^D = -0.1, [Fe/H],3y^^TD = 
—0.7, and [Fe/Hj^^gy^^ jj = —1.4, as well as 

with cr[Fe/H]D = 0.2, o-[Fo/h]td = O-^^, and 
"^[Fe/H] = 0.5, according to Allende Prieto 
et al. ""(2004, 2006). Thus, the probability 
that a star with the Galactic space veloc- 
ity components of this system (Mirabel et al. 
2001; U = -105 ± 16 km s'^ ^ = -98 ± 16 
kms~^ , W = —21 it 10 kms~^ ) and metal- 
licity [Fe/H] = 0.18 belongs to the Galactic 
halo is less than 0.1%. Moreover, the kine- 
matics alone suggest thick-disk rather than 
halo membership, although the high metal- 
licity of the secondary star favors thin-disk 
membership. However, the system could also 
have originated in a satellite galaxy. In par- 
ticular, its galactocentric orbit (Mirabel et al. 
2001) is marginally consistent with the equa- 
torial orbit of the stream of the dwarf galaxy 
in Canis Major, although its present Galac- 
tic latitude differs by more than 40° from the 
I — b distribution of the remnant of this dwarf 
galaxy, accreted by the Milky Way (Martin et 
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Fig. 8. — Chemical composition of the secondary atmosphere contaminated with nucleosynthetic products of a 
40 M0 spherically symmetric core-collapse SN explosion model (Mhb 15-16 M0) for different metallicities and 
mass-cuts, Mcut- 



al. 2004). 

The system could also have originated as 
a consequence of an encounter of an ancient 
black hole of the Galactic halo with a single 
star or a binary of two solar-type stars of the 
Galactic disk. However, this possibility is very 
unlikely due to the extremely low density of 
stars in the disk (~0.006 stars/pc^; Mihalas 
& Binney 1981). The orbit of the system in- 
tegrated backward in time never crossed the 
Galactic plane through the inner 2 kpc (Gua- 
landris et al. 2005), so high-density regions 
near the Galactic center are discarded. Porte- 
gies Zwart et al. (1997b) have modeled the 
encounter of black holes in high-density sys- 
tems. In high-density systems, where the den- 



sity of stars is ~ 4 X 10^ stars/pc^ (10^ times 
higher than in the Galactic disk) , a black hole 
spends 1.5 Gyr before it suffers a tidal cap- 
ture by a main-sequence star. In addition, the 
cases where a black hole can capture an iso- 
lated star or one star of a binary require very 
stringent constraints on the closest approach 
and impact velocity (Bcnz & Hills 1992; Hills 
1991). All these reasons make this possibility 
very unlikely. 

Therefore, in conclusion, the present loca- 
tion, velocity, and metallicity of the secondary 
star in XTE J1118-F48 suggest that the black 
hole formed in a supernova/hypernova explo- 
sion that occurred within the binary system. 
This explosive event must have either pro- 
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vided a kick to the system if it was formed 
in the thin disk, or enriched significantly the 
atmosphere of the secondary star if the system 
formed in the thick disk or halo. 

The present orbital separation between 
the compact object and the secondary star 
has been estimated to be Oc « 2.67 Rq 
(Gelino et al. 2006). Thus, the secondary 
star could have captured a significant amount 
of the ejected matter in the SN explosion that 
formed the compact object. The chemical 
composition of the secondary star may pro- 
vide information on the chemical composition 
of the progenitor of the compact object, and 
therefore on the formation region (thin disk, 
thick disk, or halo) of the binary system. We 
will now discuss the possibility that the SN 
explosion of the massive progenitor enriched 
the secondary star from different initial metal- 
licities. 

5.1. Spherical Explosion 

Gelino et al. (2006) derived a current black 
hole mass of Mbh,/ = 8.53 it 0.60 Mq and sec- 
ondary mass of M2,/ = 0.37 ±0.03 Mq. Using 
near-UV spectroscopic observations of the ac- 
cretion disk, Haswell et al. (2002) suggested 
that the material accreted onto the compact 
object is substantially CNO processed, indi- 
cating that the initial mass of the secondary 
star could have been as high as ~ 1.5 Mq. 
Hereafter we will adopt an initial secondary 
mass of M2,j = 1 Mq and a black hole mass 

of MBH,i = 8 Mq. 

A binary system such as XTE Jl 118-1-480 
will survive a spherical SN explosion if the 
ejected mass AM = Mne - Mbh.i < {Mue + 
M2,i)/2 (Hills 1983). This implies a mass 
of the He core before the SN explosion of 
Mfje < 17 Mq. Using the expressions given 
by Portegies Zwart et al. (1997a, and refer- 
ences therein), we inferred a He core radius 
of i?He ~ 2 — 3 Rq for He core masses in 
the range Mne 8.5 - 17 Mq. We will as- 



sume that the post-SN orbital separation af- 
ter tidal circularization of the orbit is in the 
range ac,i 4 — 6 Rq, since the secondary 
must have experienced mass and angular mo- 
mentum losses during the binary evolution 
until reaching its present configuration, with 
acj ~ 2.67 Rq. 

Assuming a pre-SN circular orbit and an 
instantaneous spherically symmetric ejection 
(that is, shorter than the orbital period), 
one can estimate the pre-SN orbital separa- 
tion, ao, using the relations given by van 
den Heuvel &: Habets (1984): ao = OciM/) 
where /// = (MeH.i + M2,i)/{MHe + Ma,^). 
We find oq ~ 3-5, essentially depending on 
the adopted values of Mhb and ac,i- At the 
time of the SN explosion (~5-6 Myr; Brun- 
ish &; Truran 1982), a 1 Mq secondary star, 
still in its pre-main-sequence evolution, has a 
radius i?2,i « 1.3 Rq and a convective zone 
of mass Mcz 0.652 Mq (D'Antona &; Mazz- 
itelli 1994). Thus, the amount of mass de- 
posited on the secondary can be estimated 
as mcap = AM(7ri?^^i/47rao2)/cap Mq, where 
/cap is the fraction of mass, ejected within the 
solid angle subtended by the secondary star, 
that is eventually captured. We assume that 
the captured mass, mcap, is efficiently mixed 
with the mass of the convective zone, Mcz- 

We compute the expected abundances in 
the atmosphere of the secondary star after 
the pollution from the progenitor of the com- 
pact object as in Gonzalez Hernandez et al. 
(2004). We used 40 Mq spherically symmet- 
ric core-collapse explosion models (Mfjc ~ 
15.1-16.1 Mq) for different metallicities (Z = 
0, 0.001, 0.004, 0.02) and explosion energies 
(Umeda & Nomoto, 2002, 2005; Tominaga, 
Umeda & Nomoto 2007). These models im- 
ply AM « 7-8 Mq and need small capture 
efficiencies of /cap ^ 0.1 (i.e., 10%) to in- 
crease significantly the metal content of the 
secondary star. On the other hand, the use 
of 30 Mq models (Mrb « 8.5-11.2 Mq) would 
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Fig. 9. — Left panel: Observed abundances (filled circles with error bars) in comparison with the expected abun- 
dances in the secondary star after having captured 17.5% of the matter ejected within the solid angle subtended by the 
secondary from a metal-poor {Z = 0.001) 40 M0 spherically symmetric supernova explosion (Mho = 15.8 Mq) with 
Ek = 10""^ erg for two different mass-cuts, Mcut = 2.51 Mq (solid line with open circles) and Mcut = 3.01 Mq (dashed- 
dotted line with open circles). The initial abundances of the secondary star were adopted for the average abundances 
of halo stars with [Fe/H] = —1.4 ± 0.2, and the initial orbital distance was Oc.i ~ 6 R©. Right panel: same as 
left panel but for a spherically symmetric hypernova explosion {Ek = 30 x 10^^ erg) for two different mass-cuts, 
Mcut ~ 4.05 M0 (solid line with open circles) and Af^ut = 5.03 Mq (dashed-dotted line with open circles). 



require /cap ~ 0.9-1. These models would also 
provide a different mass fraction of each ele- 
ment at each value of the mass-cut (the mass 
that initially collapsed forming the compact 
remnant). For more details in the models, see 
Tominaga, Umeda & Nomoto (2007). 

The explosion energy, Ek = 1 x 10^^ erg 
and Ek = (20 - 30) x 10^^ erg for the su- 
pernova (SN) and hypernova (HN) models 
(respectively), is deposited instantaneously in 
the central region of the progenitor core to 
generate a strong shock wave. The subse- 
quent propagation of the shock wave is fol- 
lowed through a hydrodynamic code (Umeda 
&: Nomoto 2002, and references therein). In 
our simple model, we have assumed different 
mass-cuts, fallback masses, and a mixing fac- 
tor of 1 which assumes that all fallback matter 
is well mixed with the ejecta. The amount of 
fallback, Mfaii, is the difference between the fi- 
nal remnant mass, M^n^i, and the initial rem- 
nant mass of the explosion, Mcut- We should 
recall here the ejected mass, AM, which is 



equal to Mhc — M^u^i, where Mne is the mass 
of the He core. 

We use SN/HN models to provide us with 
the yields of the explosion before radiative 
decay of element species. We then compute 
the integrated, decayed yields of the ejecta by 
adopting a mass-cut and by mixing all the ma- 
terial above the mass-cut. Finally, we calcu- 
late the composition of the matter captured 
by the secondary star, and we mix it with the 
material of its convective envelope. 

In Fig. [8] we show the expected abun- 
dances of the secondary star after contami- 
nation from the nucleosynthetic products of 
the SN explosion {Ek = 1 x 10^^ erg) of 
Mne ~ 15-16 Mq progenitor stars. The ini- 
tial abundances of the secondary star have 
been estimated from the average abundances 
of halo stars with a metallicity of [Fe/H] 
« -2.2 (from Cayrel et al. 2004) and [Fe/H] 

-1.4 (from Jonsell et al. 2005), thick-disk 
stars with [Fe/H] ~ —0.7 (from Jonsell et al. 
2005), and thin-disk stars with [Fe/H] « 
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Fig. 10. — Left panel: Observed abundances (filled circles with error bars) in comparison with the expected 
abundances in the secondary star after having captured the 15.5% of the matter ejected within the solid angle 
subtended by the secondary from a metal-poor {Z — 0.004) 40 Mq spherically symmetric supernova explosion (Mhc ~ 
15.8 Mq) with Ek = lO^'^ erg for two different mass-cuts, Mcut ~ 2.50 M© (solid line with open circles) and 
A/cut = 3.03 Mq (dashed-dotted line with open circles). The initial abundances of the secondary star were adopted 
for the average abundances of thick-disk stars with [Fe/H] = — 0.7±0.2, and the initial orbital distance was ac,i « 6 Kq. 
Right panel: same as left panel but for a spherically symmetric hypernova explosion {Ek = 30 x 10^^ erg) for two 
different mass-cuts, Mcut ~ 3.03 Mq (sohd line with open circles) and Aleut = 4.02 Mq (dashed-dotted line with open 
circles) . 



and [Fe/H] ^ 0.18 (from Gilli et al. 2006). 
For each simulation at each metahicity, we 
have fixed the /cap factor at a value that al- 
lows us to approximately match the observed 
aluminum abundance. 

We should remark that for a given model at 
a given metallicity, once the capture efficiency 
is fixed, the aluminum (Al) abundance in the 
secondary star hardly depends on the mass- 
cut. Thus, the fact that the Al abundance 
in Fig. [8] is almost constant with metallicity is 
because we have adopted different /cap factors 
for different metallicities. These /cap factors 
are higher for lower metallicities since greater 
amounts of captured material are needed to 
achieve the observed Al abundance in the sec- 
ondary star. For other figures in this paper, 
the /cap factor was changed until an abun- 
dance of [Al/H]« 0.5 dex was obtained, com- 
patible with the observed abundance within 
the uncertainties. 

From Fig. [8] one can see that while Mg and 



Al remain significantly enhanced (above solar 
abundances) for all mass-cuts at all metallici- 
ties, Si, Ca, Ti, Fe, and Ni are quite sensitive 
to the mass-cut of the model. Thus, for metal- 
licities below —1.4 dex, only Mcut ^ 3 Mq is 
able to enhance sufficiently the Ti, Fe, and Ni, 
whereas Ca and Si remain quite overabundant 
for Mcut ^ 4 Mq and 5 Mq (respectively) 
at all metallicities. We should note that 
for metallicities [Fe/H] < —1.4 dex (i.e., 
Z < 0.001), the expected abundances in the 
secondary star do not depend on the initial 
abundances but on the SN yields. Thus, we 
found different expected abundances of the 
secondary between [Fe/H] ~ —1.4 dex and 
[Fe/H] « —2.2 dex because we used the SN 
Z = 0.001 model and the SN Z = model, 
respectively. 

5.2. Formation in the Halo 

The kinematics of the system, at least U 
and V, resemble those of halo stars, and there- 
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Table 5 

Metal-Poor Supernova/Hypernova Explosion Models in XTE J1118-I-480 



[X/H] EXPECTED'^ 



Supernova Hypernova 



element 


[X/H] OBSERVED'^ [X/HJo^ 


-^cut.low 


" Mcut.up 


A^cut,low 










40 Mq explosion model of Z = 


= 0.001 


Mg 


0.35 


-1.08 


0.83 


0.85 


0.75 


0.79 


Al 


0.60 


-1.27 


0.52 


0.53 


0.39 


0.42 


Si 


0.37 


-1.13 


0.88 


0.89 


1.08 


1.09 


Ca 


0.15 


-1.12 


0.74 


0.75 


1.09 


1.01 


Ti 


0.32 


-1.15 


0.14 


0.11 


0.54 


0.12 


Fe 


0.18 


-1.46 


0.47 


0.14 


0.63 


-0.22 


Ni 


0.30 


-1.48 


1.09 


-o.r52 


0.14 


-0.81 


o 




-0.76 


0.75 


0.77 


0.69 


0.73 


c 




-0.86 


-0.54 


-0.53 


-0.60 


-0.58 








40 M0 


explosion model with Z 


= 0.004 


Mg 


0.35 


-0.41 


0.58 


0.59 


0.61 


0.64 


Al 


0.60 


-0.46 


0.51 


0.52 


0.41 


0.44 


Si 


0.37 


-0.50 


0.62 


0.60 


0.90 


0.93 


Ca 


0.15 


-0.54 


0.57 


0.42 


0.86 


0.87 


Ti 


0.32 


-0.50 


0.16 


-0.18 


0.56 


0.27 


Fe 


0.18 


-0.75 


0.26 


-0.52 


0.67 


0.07 


Ni 


0.30 


-0.73 


0.08 


-0.53 


0.43 


-0.38 


O 




-0.26 


0.71 


0.73 


0.67 


0.69 


c 




-0.35 


0.10 


0.11 


0.06 


0.08 








30 M© 


explosion model with Z 


= 0.004 


Mg 


0.35 


-0.41 


0.60 


0.62 


0.57 


0.62 


Al 


0.60 


-0.46 


0.52 


0.54 


0.40 


0.45 


Si 


0.37 


-0.50 


0.77 


0.54 


1.02 


0.95 


Ca 


0.15 


-0.54 


0.72 


0.17 


1.03 


0.82 


Ti 


0.32 


-0.50 


0.18 


-0.41 


0.66 


0.01 


Fe 


0.18 


-0.75 


0.14 


-0.71 


0.71 


-0.32 


Ni 


0.30 


-0.73 


-0.43 


-0.54 


0.40 


-0.40 


O 




-0.26 


0.75 


0.78 


0.69 


0.74 


c 




-0.35 


0.02 


0.03 


-0.04 


-0.01 



'^Observed abundances of the secondary star in XTE ,111184-480. 

'^Initial abundances assumed for the secondary stax in XTE J1118-I-480, see text. The initial 
C and O abundances of the thick-disk model ([X/H]q = —0.7) were adopted from Ecuvillon 

et al. (2004, 2006). 

"^A^cutjiow a^nd A^cut.up are the lower and upper mass-cut adopted in the model computation. 
See the exact value in the captions of Figs. 9-11. 

'^Expected abundances of the secondary star. 

Note. — Expected abundances in the secondary atmosphere contaminated with nucleosyn- 
thetic products of metal-poor explosion models for two different explosion energies and mass- 
cuts, presented in Figs. 9-11. 
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fore a significant kick during the black hole 
formation process appears unnecessary. Thus, 
a spherically symmetric SN explosion would 
provide the desired kick velocity to match 
the current velocity components of the sys- 
tem from initial velocities similar to those of 
halo stars. However, due to the extremely low 
metallicities of halo stars, it is required that 
the secondary captured enough matter from 
the ejecta to reach the current abundances. 

In Fig. [9] we compare the observed abun- 
dances with the expected abundances in the 
secondary star for two different explosion 
energies and mass-cuts using a metal-poor 
40 Mq progenitor model. In these simula- 
tions, the mass-cut is sampled according to 
the mass bins given in the explosion models, 
from ~ 1.5 Mq to 8 Mq in steps of ~ 0.5 and 
1 Mq . The adopted mass-cuts shown in Fig. [9] 
were selected from those which provided bet- 
ter fits to the observed abundances. The ini- 
tial abundances of the secondary have been es- 
timated from the average abundances of halo 
stars with [Fe/H] ~ —1.4 from Jonsell et al. 
(2005). The value of /cap has been adjusted 
until the observed aluminum abundance was 
roughly reproduced, providing /cap ~ 0.17. 
The parameters of the explosion models used 
in Fig. [9] - [13] are given in Table [H In Ta- 
ble [5] we show the expected abundances of 
the secondary star after contamination from 
metal-poor SN/HN explosion models. 

With a supernova model, shown in the left 
panel of Fig. [9l it is not possible to recover 
the observed abundances because these mod- 
els produce too much Mg and Ca regardless of 
the Fe abundance obtained; in addition, the 
Ni abundance is strongly dependent on the 
mass-cut. The hypernova model, displayed in 
the right panel, makes it even more difficult to 
fit the observed abundances since this model 
creates too much Ca relative to Fe and Ni. 
As inferred from Fig. [51 a model with initial 
abundances at [Fe/H] ^ —2.2 (Cayrel et al. 



2004) would also not be successful in repro- 
ducing the observed abundances due to the 
tremendous and different sensitivity of each 
element to the mass-cut. 

In conclusion, neither of these very metal- 
poor models is able to reproduce the ob- 
served abundances. We also try to fit the 
observed abundances using a metal-poor 
30 Mq explosion model, but the agreement is 
even worse than for 40 Mq models, again due 
to the strong sensitivity of the each element 
abundance to the mass-cut. The comparison 
of the observed abundances with SN yields 
allows us to rule out a Galactic halo origin for 
this black hole binary. 

5.3. Formation in the Thick Disk 

The space-velocity components of the sys- 
tem are comparable to those of thick-disk 
stars, and its present location 1.6 kpc above 
the Galactic plane is slightly higher that the 
scale height of thick-disk stars (~ 0.8-1.3 kpc; 
Reyle & Robin 2001; Chen 1997); thus, a 
strong kick during the SN explosion would not 
be required. A spherically symmetric SN ex- 
plosion of a 15.8 Mq and 11 Mq He core pro- 
vides an impulse of ~ 60 and ~ 20 km s~^ , re- 
spectively. However, an enrichment from the 
typical abundances of thick-disk stars would 
have been necessary. In the simulations, the 
initial abundances were assumed to be the av- 
erage values of thick-disk stars with [Fe/H] 
^ -0.7 from Jonseh et al. (2005). 

In Fig. [To] we compare the expected abun- 
dances from a 40 Mq explosion model for two 
energies and mass-cuts (see also Table [5]). 
The left panel shows the expected abundances 
from a SN model which seems to better ap- 
proach the observed abundances than those 
of halo-like metallicities. As in the previ- 
ous figure, the /cap factor was changed un- 
til an abundance of [Al/H]~ 0.5 dex was ob- 
tained, compatible with the observed abun- 
dance within the uncertainties. For the mass- 
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Table 6 

Metal-Rich Supernova/Hypernova Explosion Models in XTE J1118+480 



[X/H] EXPECTEpd 

ELEMENT [X/H] OBSERVED'' [X/H]^^ M^ut.iow'^ Mcut.up Meut,io„ Mcut,, 

Si)li(>ri<.*al cxpiusioii uio<[(4 of Z = 0,02 
Supernova Hypernova 



Mg 


0.35 


0.17 


0.32 


0.30 


0.32 


0.30 


Al 


0.60 


0.29 


0.50 


0.50 


0.48 


0.50 


Si 


0.37 


0.12 


0.22 


0.14 


0.28 


0.14 


Ca 


0.15 


0.02 


0.09 


0.02 


0.15 


0.02 


Ti 


0.32 


0.20 


0.24 


0.21 


0.34 


0.21 


Fe 


0.18 


0.18 


0.30 


0.18 


0.32 


0.18 


Ni 


0.30 


0.13 


0.35 


0.15 


0.32 


0.15 


O 




0.18 


0.33 


0.33 


0.33 


0.33 


c 




0.11 


0.16 


0.16 


0.15 


0.16 



Aspherical explosion model with Z = 0.02 
Angle" = 0° - 15° Angle = 0° - 90° 



Mg 


0.35 


0.17 


0.34 


0.32 


0.34 


0.35 


Al 


0.60 


0.29 


0.51 


0.48 


0.50 


0.51 


Si 


0.37 


0.12 


0.24 


0.18 


0.35 


0.26 


Ca 


0.15 


0.02 


0.05 


0.02 


0.22 


0.11 


Ti 


0.32 


0.20 


0.21 


0.20 


0.41 


0.35 


Fe 


0.18 


0.18 


0.18 


0.18 


0.28 


0.22 


Ni 


0.30 


0.13 


0.13 


0.13 


0.34 


0.23 


O 




0.18 


0.37 


0.35 


0.38 


0.39 


c 




0.11 


0.12 


0.13 


0.12 


0.13 



''Observed abundances of the secondary star in XTE J1118+480. 

*^ Initial abundances assumed for the secondary star in XTE J1118+480, see text. The initial 
C and O abundances of the metal-rich models were adopted from Ecuvillon et al. (2004, 2006). 

°Afcut,iow and Mcut, up are the lower and upper mass-cut adopted in the model computation. 
See the exact value in the captions of Figs. 12 and 13. 

•^Expected abundances of the secondary star. 

''Angular range, measured from the equatorial plane, in which all the ejected material in 
the explosion has been completely mixed for each velocity point. 

Note. — Expected abundances in the secondary atmosphere contaminated with nucleosyn- 
thctic products of metal-rich explosion models for two different mass-cuts and symmetries, 
presented in Figs. 12 and 13. 
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Fig. 11. — Left panel: Observed abundances (filled circles with error bars) in comparison with the expected 
abundances in the secondary star after having captured 90% of the matter ejected within the solid angle subtended by 
the secondary from a metal-poor {Z = 0.004) 30 Mq spherically symmetric supernova explosion (Mhc = 11 Mq) with 
Ek = 10""^ erg for two different mass-cuts, Afcut = 2.52 Mq (solid line with open circles) and Mcut = 3.04 Mq (dashed- 
dotted line with open circles). The initial abundances of the secondary star were adopted for the average abundances 
of the thick-disk stars with [Fe/H] = —0.7 ± 0.2, and the initial orbital distance was ac,i ~ 6 R©. Right panel: same 
as left panel but for a spherically symmetric hypernova explosion {Ek = 20 x 10^^ erg) for two different mass-cuts, 
Mcut = 3.04 Mq (solid line with open circles) and M^ut = 4.15 Mq (dashed-dotted line with open circles). 



cut of 2.50 Mq, the model roughly fits the ob- 
served abundances except for Ca, which ap- 
pears enhanced in the model by a factor of 
2.6. Better agreement with Ca could be ob- 
tained if the /cap factor is lowered, but then 
Al, Ni, and Ti would not match their observed 
abundances. If we increase the mass-cut up 
to 3.03 Mq, the Ca abundance only decreases 
0.15 dex whereas Ti, Fe, and Ni decrease by 
a factor of 2, 6, and 4 (respectively), which 
makes the model unable to fit the observed 
abundances. The hypernova case offers worse 
results since Si and Ca are no longer repro- 
duced, and Fe and Ni cannot be fitted at the 
same time. 

For this metallicity we have also explored 
the 30 Mq explosion model. The results are 
displayed in Fig. [TTl In the left panel the SN 
model is compared with the observations. For 
the lower mass-cut, at 2.52 Mq, the model 
provides too high Si and Ca abundances and 
too low Ni abundance, whereas for a mass- 
cut of 3.04 Mq, the observed Ti, Fe, and Ni 



abundances are too high in comparison with 
the model predictions. In the hypernova case 
(right panel of Fig. [TT|) . at these low mass- 
cuts the model produces too large amounts of 
Si and Ca; in addition, Ti, Fe, and Ni are too 
sensitive to the location of the mass-cut. 

Despite the fact that none of the explosion 
models explored is able to fairly reproduce 
the observed abundances in the secondary 
star, all of the models require vigorous mix- 
ing between the fallback matter and the fi- 
nal ejected matter (Kifonidis et al. 2000). For 
instance, in the 15.8 Mq He core model, all 
of the material above the mass-cut placed at 
2.50 Mq should be well mixed in order to con- 
vey heavy elements like Fe and Ni to the outer 
layers of the explosion, which might make 
more unlikely a Galactic thick-disk origin for 
XTE J1118-F480. 

5.4. Formation in the Thin Disk 

In this scenario, the system must have ac- 
quired an impulse during the formation of the 
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Fig. 12. — Left panel: Observed abundances (filled circles witii error bars) in comparison with the expected abun- 
dances in the secondary star after having captured the 5% of the matter ejected within the solid angle subtended by 
the secondary from a metal-poor {Z — 0.02) 40 M0 spherically symmetric supernova explosion (Mho = 15.1 M0) with 
Ek ~ 10^^ erg for two different mass-cuts, Mcut = 1.46 Mq (solid line with open circles) and Mcut = 8.09 Mq (dashed- 
dotted line with open circles) . The initial abundances of the secondary star were adopted for the average abundances 
of thin-disk stars with [Fe/H] = 0.0 ± 0.2. Right panel: same as left panel but for a spherically symmetric hyper- 
nova explosion {Ek = 30 x 10^^ erg) for two different mass-cuts, Mcut ~ 1-74 Mq (solid line with open circles) and 
Mcut = 8.09 M0 (dashed-dotted line with open circles). 



black hole that pushes it up from a Galactic 
plane orbit to the current halo orbit. The sys- 
tem should have been accelerated to a peculiar 
space velocity of ~ 180 km s~^ (Gualandris et 
al. 2005) to reach its present location, requir- 
ing an asymmetric kick. It has been suggested 
that such kicks can be imparted during the 
birth of nascent neutron stars, due to asym- 
metric mass ejection and/or an asymmetry in 
the neutrino emission (Lai et al. 2001, and 
references therein). 

Podsiadlowski et al. (2002) proposed that 
the black hole in Nova Sco 1994 could have 
formed in a two-stage process where the ini- 
tial collapse led to the formation of a neu- 
tron star accompanied by a substantial kick 
and the final mass of the compact remnant 
was achieved by matter that fell back af- 
ter the initial collapse. However, the black 
hole mass in Nova Sco 1994 is estimated to be 
~ 5.4 M0 (Beer & Podsiadlowski 2002), while 
the black hole in XTE J1118-F480 has a mass 
of ~ 8 M0 which would require a fallback 



mass of ^ 6.6 Mq if we assume ~ 1.4 Mq for 
a canonical neutron star. This scenario might 
take place in the context of collapsar mod- 
els where the black hole would be formed in 
a mild explosion and substantial fallback (up 
to ~ 5 Mq is expected) as proposed by Mac- 
Fadyen et al. (2001). Asymmetric mass ejec- 
tion would relax this requirement, providing 
enough impulse to the system to be launched 
into its present orbit from the Galactic plane. 

In Fig. [T2t we compare the expected abun- 
dances from the explosion of a 15.1 Mq He 
core with the observed abundances of the sec- 
ondary star (see also Table [6]). As in the pre- 
vious figures, the /cap factor was changed un- 
til an abundance of [Al/H]~ 0.5 dex was ob- 
tained, compatible with the observed abun- 
dance within the uncertainties. The initial 
abundances were assumed to match the av- 
erage values of disk stars with similar iron 
content, which are provided in Table [3l For 
both the SN model (left panel) and the hy- 
pernova model (right panel), the observed el- 
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Fig. 13. — Left panel: Observed abundances (filled circles with error bars) in comparison with the expected 
abundances in the secondary star after having captured the 5% of the matter ejected within the solid angle subtended 
by the secondary from a non-spherically symmetric supernova explosion of Ek ~ W^^ erg for two different mass-cuts, 
Mcut = 2.41 Mq (solid line with open circles) and Mcut = 7.56 Mq (dashed-dotted line with open circles). This 
model corresponds to the matter ejected in the equatorial plane of the primary where we assumed that the secondary 
star is located (for more details, see Gonzalez Hernandez et al. 2005b). Right panel: same as left panel but in 
this model we have assumed complete lateral mixing, where all the material within given velocity bins is completely 
mixed. Two simulations are shown for two different mass-cuts, Mcut ~ 2.41 Mq (solid line with open circles) and 
Mcut = 7.56 M© (dashed-dotted line with open circles). 



ement abundances can be reproduced for all 
mass-cuts. This means that neither substan- 
tial fallback nor an efficient mixing process is 
needed. For mass-cuts above ~ 3 M©, very lit- 
tle Si, Ca, Ti, Fe, and Ni is ejected; thus, the 
expected abundances of the model essentially 
reflect the initial abundances of the secondary 
star. In contrast, Mg and Al are not sensitive 
to the mass-cut and are slightly enhanced due 
to the capture of enriched material in the SN 
explosion. 

We also investigate a model with solar ini- 
tial abundances for the secondary star, i.e., 
[X/H]q = 0, and the same explosion model of 
solar metallicity {Z = 0.02), which would fit 
all the element abundances if the mass-cut is 
low enough (Mcut ^ 3 Mq) and mixing is so 
efficient that significant amounts of elements 
which form in the inner layers of the explo- 
sion (such as Ti, Ni, and Fe) are present in 
the ejecta. 



5.5. Non-Spherical Explosion: Forma- 
tion in the Thin Disk 

In this section, the thin-disk scenario is 
studied using non-spherically symmetric ex- 
plosion models from Maeda et al. (2002). 
The chemical composition of the ejecta in 
a non-spherically symmetric SN explosion is 
strongly dependent on direction. In particu- 
lar, if we assume that the jet is collimated per- 
pendicular to the orbital plane of the binary, 
where the secondary star is located, elements 
such as Ti, Fe, and Ni are mainly ejected in 
the jet direction, while O, Mg, Al, Si, and S 
are preferentially ejected near the equatorial 
plane of the helium star (Maeda et al. 2002). 

In Fig. [13] we compare the predicted abun- 
dances in the secondary star after pollution 
from an aspherical explosion model of a metal- 
rich progenitor having a 16 Mq He core (see 
also Table [6]) . The initial abundances of the 
secondary were extracted from the average 
values of stars of the solar neighborhood with 
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similar iron content (see Table [3|). The left 
panel reflects the composition of the mate- 
rial ejected in the equatorial plane while in 
the right panel we have considered complete 
lateral mixing (Podsiadlowski et al. 2002) — 
that is, the ejected matter is completely mixed 
within each velocity bin. This mixing pro- 
cess is done with the decayed yields of the 
model, and after that, all the material above 
the mass-cut is mixed and used to calculate 
the composition of the ejected matter. The 
observed abundances are better reproduced 
if complete lateral mixing is considered since 
this process tends to enhance all of the ele- 
ment abundances at all mass-cuts. However, 
the equatorial model (left panel) also pro- 
vides good fits to the observed abundances. 
It should be noted that in this model, only 
the material ejected in the equatorial plane 
is captured and therefore, only Mg, Al, and 
Si are significantly enhanced with respect to 
the initial abundances. A model with solar 
initial abundances for the secondary star, i.e. 
[X/H]q = 0, and the same explosion model 
of solar metallicity (Z = 0.02) was also in- 
spected, providing the same result except for 
the equatorial model which produces too low 
abundances of Ti, Fe, and Ni in comparison 
with the observations. 

Both aspherical explosion models with ini- 
tial abundances equal to the average values 
of thin-disk stars show little dependence on 
the mass-cut, and since the equatorial model 
can reproduce the observed abundances, ex- 
tensive mixing processes are not required. In 
addition, the non-spherically symmetric ejec- 
tion of the mass in the explosion could provide 
the kick that the system needs to change its 
orbit from the Galactic plane to the current 
orbit. Therefore, the element abundances in 
the secondary star and the kinematics of this 
system strongly suggest that the binary sys- 
tem XTE Jl 118+480 formed in the Galactic 
disk (probably the thin disk) and it was then 



kicked towards the Galactic halo, most prob- 
ably by asymmetric mass ejection in an asym- 
metric supernova/hypernova explosion that 
gave rise to the black hole in this system. 

The elements O and C, for which we pro- 
vide the expected abundances in Tables [5] 
and [6l could be studied in future investiga- 
tions, probably from CO and OH bands in 
the near-IR. This would help to recognize the 
operation of the CNO cycle on the surface 
the secondary star, proposed by Haswell et 
al. (2002), and possible processes of rotation- 
induced mixing during the evolution of the 
massive star. 

6. Conclusions 

We have presented Keck II/ESI medium- 
resolution spectroscopy of the black hole bi- 
nary XTE J1118+480. The individual spec- 
tra of the system allowed us to derive an or- 
bital period of P = 0.16995 ± 0.00012 d and 
a radial velocity semiamplitude of the sec- 
ondary star of K2 = 708.8 it 1.4 kms^"*^. The 
implied updated mass function is /(M) = 
6.27 lb 0.04 Mq , consistent with (but more 
precise than) previous values reported in the 
literature. Inspection of the high-quality av- 
eraged spectrum of the secondary star pro- 
vides a rotational velocity of v sin i = lOOt^i 
kms~^, and hence a binary mass ratio q = 
0.027 ± 0.009. The derived radial velocity, 
7 = 2.7 lb 1.1 kms^"*^, of the center of mass 
of the system agrees, at the 3cr level, with the 
results of previous studies. 

We have performed a detailed chemical 
analysis of the secondary star. We applied a 
technique that provides a determination of the 
stellar parameters, taking into account any 
possible veiling from the accretion disk. We 
find Teff = 4700 + 100 K, log[c//cm s^] = 4.6 ± 
0.3, [Fe/H] = 0.18 ib 0.17, and a disk veiling 
(defined as Fdisk/i^totai) of ~40% at 5000 A, 
decreasing toward longer wavelengths. 



22 



We have provided further details on the 
abundances of Mg, Al, Ca, Fc, Ni, and Li 
already reported by Gonzalez Hernandez et 
al. (2006), and we determined new element 
abundances of Si and Ti. The chemical abun- 
dances are typically higher than solar, and in 
some cases they are slightly enhanced (e.g., 
Mg, Al, and Si) in comparison with the abun- 
dances of these elements in stars of the solar 
neighborhood having similar iron content. 

The present location and kinematics of this 
binary system had suggested that it could 
have originated in the Galactic halo. How- 
ever, the chemical abundances strongly indi- 
cate that the black hole formed as a conse- 
quence of a supernova/hypernova explosion 
that occurred within the binary system. This 
explosive event must have either provided a 
kick to the system if it was formed in the thin 
disk or enriched significantly the atmosphere 
of the secondary star if the system formed in 
the thick disk or halo. 

We have explored a variety of super- 
nova/hypernova explosion models for differ- 
ent metallicities, He core masses, and ge- 
ometries. We compared the expected abun- 
dances in the secondary star after contam- 
ination from nucleosynthetic products from 
different initial metallicities of the secondary 
star (-2.2 < [Fe/H] < 0.2), to investigate the 
formation region in the Galactic halo, thick 
disk, or thin disk. Metal-poor explosion mod- 
els {Z = OandZ = 0.001) were not able to fit 
the observed abundances since they produce 
inappropriate ratios between a-elements and 
iron-peak elements, and they are extremely 
sensitive to the adopted mass-cut. This com- 
parison probably rules out an origin in the 
Galactic halo for this black hole binary. 

For the thick-disk scenario, we carefully in- 
spected the model predictions, and although 
they provide better fits to the observed abun- 
dances, they require substantial fallback (up 
to 5.5 Mq ) and very efficient mixing processes 



between the inner layer of the explosion and 
the ejecta. We thus conclude that this sce- 
nario is unlikely. 

Metal-rich spherically symmetric models 
for the thin-disk scenario were able to fairly 
reproduce the observed abundances, although 
they do not easily provide the energy required 
to launch the system from the Galactic plane 
to its current halo orbit. 

Finally, non-spherically symmetric mod- 
els produce excellent agreement with the ob- 
served element abundances in the secondary 
star without invoking extensive fallback and 
mixing. In addition, asymmetric mass ejec- 
tion would naturally provide the kick to ex- 
pel this binary system from its birth place in 
the Galactic thin disk, which seems to be the 
most plausible explanation for the origin of 
this halo black hole X-ray binary. 
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